Collisions between tunable halo dimers: 
exploring an elementary four-body process with identical bosons 
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We study inelastic collisions in a pure, trapped sample of Feshbach molecules made of bosonic cesium atoms 
in the quantum halo regime. We measure the relaxation rate coefficient for decay to lower-lying molecular states 
and study the dependence on scattering length and temperature. We identify a pronounced loss minimum with 
varying scattering length along with a further suppression of loss with decreasing temperature. Our observations 
provide insight into the physics of a few-body quantum system that consists of four identical bosons at large 
values of the two-body scattering length. 

PACS numbers: 21.45.-v, 33.15.-e, 34.50.Cx, 37.10.Pq 
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Few-body quantum systems in halo states exhibit unique 
properties 1 1 ] . A quantum halo is a very weakly bound state 
whose wave function extends far into the classically forbidden 
range. Halo systems are much larger than one would expect 
from the characteristic interaction range of their constituents. 
Many examples for halo states are known in nuclear physics, 
with the deuteron being a prominent example J2]. In molec- 
ular physics, the He dimer has for many years served as the 
prime example of a halo state 13,0]. 

A particular motivation to study halo states is given by the 
concept of universality in few-body systems |5|]. Since, in 
a halo state, short-range details of the interaction become ir- 
relevant, the system is described by very few parameters and 
shows universal behavior in its low-energy observables. A 
halo dimer is the elementary two-body halo system. Here the 
only relevant length scale is given by the scattering length a, 
which describes the s-wave interaction between its two con- 
stituents. The size of the halo dimer is directly related to a 
and the binding energy is Eb = ft 2 / (2\ia 2 ), where jl is the re- 
duced mass. For three-body halo systems, universal Efimov 
states can exist |0, 0]. Here one additional parameter is re- 
quired to fully describe the system in the universal limit. A 
natural further step is to investigate the universal physics of 
systems composed of four identical bosons. The fundamen- 
tal properties of such systems are unexplored terrain, with the 
existence of an additional four-body parameter [8, 9], the scal- 
ing and threshold behavior, and the binding energies of four 
body-states [9] being open issues. 

In the field of ultracold gases, the Feshbach association 
technique H 1 Oil has provided experimentalists with unprece- 
dented possibilities to create and study halo dimers. In the 
case of certain Feshbach resonances, a considerable range of 
universality exists, where halo dimers can be conveniently 
controlled by a magnetic bias field to vary their binding en- 
ergy and size. Such tunable halo dimers are unique probes to 
explore quantum phenomena related to universality. A binary 
collision between two halo dimers can be seen as an elemen- 
tary four-body process. For a special kind of halo dimers such 
a four-body process has already attracted considerable atten- 
tion: halo dimers made of fermionic atoms in different spin 



states allow to create molecular Bose-Einstein condensates 
and to study the crossover to a fermionic superfluid [11]. Here 
a key point is the Pauli suppression effect that, in combination 
with the halo nature of the dimer, leads to stability against de- 
cay into lower-lying molecular states and favors elastic pro- 
cesses fl2ll . 

In this Letter, we study binary collisions in a pure, trapped 
sample of tunable halo dimers made of bosonic atoms. Halo 
dimers of this class have so far received much less experimen- 
tal attention than their fermionic counterparts, although they 
represent an important link to universal few-body phenom- 
ena in systems of few interacting bosons; for three particles, 
an early example is the predicted atom-dimer "Efimov" reso- 
nances W13N . Halo dimers of bosonic atoms have been realized 
in ultracold gases of 85 Rb and 133 Cs 0, 03] and properties 
of the individual dimers, like binding energies, magnetic mo- 
ments, and spontaneous decay rates, have been measured. In 
contrast, their collision properties have remained unexplored 
terrain. Because of the absence of a Pauli suppression effect, 
substantial inelastic decay to lower-lying molecular states can 
be expected. The observation of loss serves as a probe for 
dimer-dimer interactions flL fToh. 

Our experiments are performed with 133 Cs atoms, which 
represent an excellent system to study few-body physics with 
bosons at large scattering lengths [7] because of the unique 
scattering properties [16]. In the lowest spin state at low mag- 
netic fields, one finds a broad entrance-channel dominated s- 
wave Feshbach resonance [OHO] along with an extraordinar- 
ily large background scattering length. The scattering proper- 
ties of ultracold atoms are governed by the last bound s-wave 
state below the dissociation threshold as displayed in Fig.Q] In 
a wide magnetic-field range, this state carries a quantum halo 
character, where the two-body scattering length far exceeds 
the classical interaction range of the van der Waals potential, 
r vdw (~ lOO^o), and the binding energy is much smaller 
than the corresponding £ v dw (~ /z x 2.7 MHz) 118]. Here ao 
is Bohr's radius. 

Our experimental procedure to produce an optically trapped 
sample of tunable halo dimers involves several stages. We 
initially prepare ultracold trapped 133 Cs atoms in their abso- 
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FIG. 1: (color online) Weakly bound Cs2 molecular states at low 
magnetic fields, (a) Binding energies of the relevant Cs2 dimer states. 
The zero-energy level is the s-wave threshold of two colliding Cs 
atoms in their absolute hyperfine ground state sub-level, (b) s-wave 
scattering length a versus the magnetic field (see text). The shaded 
regions in (a) and (b) indicate the non-halo regime with > £ v dW> 
and a < r v( jw> respectively. 



lute hyperfine ground state sublevel |F,mir) = |3,3), similarly 
to Ref. lfl5ll . The atoms are optically trapped by two crossed 
1064-nm laser beams with waists of about 250 /im and 36 
jUm, while a magnetic field levitates the atoms against gravity 
The levitation field ensures an optimized evapora- 
tive cooling of the atoms, which is realized by lowering the 
optical power in the trapping beams. We stop the cooling just 
before the onset of Bose-Einstein condensation to avoid too 
high atomic densities. At our lowest temperature of 20 nK, 
we obtain about 1.5 x 10 5 non-condensed atoms. 

In the next stage, we create the halo dimers by Feshbach 
association. Here, the application of the levitation field is 
not appropriate since atoms and dimers in general have dif- 
ferent magnetic moments. This leads to special requirements 
for the trap design. On the one hand, we need a sufficiently 
high optical gradient in the vertical direction to hold the atoms 
and dimers against gravity. On the other hand, we want to 
avoid the high density of a tight trap, which causes fast losses 
driven by atom-dimer collisions. These two requirements can 
be simultaneously fulfilled by using an elliptic trap potential 
with weak horizontal confinement and tight confinement in 
the vertical direction. Shortly before molecule production, 
we adiabatically convert the levitated trap to a non-levitated 
trap by simultaneously changing the levitation field, the opti- 
cal power, and the trap ellipticity. The latter is modified by a 
rapid spatial oscillation of the 36-/1 m waist beam in the hor- 
izontal plane with the use of an acousto-optic modulator at a 
frequency of about 100 kHz, which greatly exceeds the typi- 
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FIG. 2: (color online) Measurements on the trapped molecular sam- 
ple, (a) Comparison between the atomic and molecular temperatures 
r a and T m at 35G, where ~ h x 87 kHz. The solid line indi- 
cates equal temperatures for atoms and molecules, (b) Number of 
halo dimers N as a function of the holding time in trap at 28.3 G 
(a = 500ao, squares) and at 45.6 G (a = 900ao, circles). The solid 
lines are fits to the data (see text). 



cal trap frequencies, and thus creates a time-averaged optical 
potential [20» |2ll |2^]. The adiabatic change of the trap shape 
is set in a way to keep the peak density, the temperature, and 
thus the phase-space density constant 02311 . Right after con- 
verting the trap, we associate the halo dimers by sweeping the 
magnetic field across the 200 mG wide d-wave Feshbach res- 
onance located at approximately 48 G 11511 . (see Fig.Q]). 

The molecular temperature can be set by adjusting the tem- 
perature of the initial atomic sample. To selectively measure 
the atomic and molecular temperatures, we spatially separate 
the two components with the Stern-Gerlach technique, and we 
perform a subsequent time-of-flight imaging, as described in 
Refs. [l9ll . In Fig.|2a), we compare the atomic and the 
molecular temperatures for a wide range of the initial tem- 
peratures. In spite of some overall heating in the conversion 
process, we observe that the dimers and the atoms have the 
same translational temperature in the trap. This observation 
may indicate an elastic interaction between atoms and dimers 
on the 10 ms time scale of our preparation sequence and it 
allows us to conveniently use the atomic sample to determine 
the temperature of the molecular gas. 

In the last step of the preparation sequence, we selectively 
remove the atoms from the dipole trap. This is done by us- 
ing a double-resonant purification scheme, which combines 
microwave excitation with resonant light I124I1 . similarly to 
Ref. II25I1 . By absorption imaging after the Stern-Gerlach tech- 
nique we verify that no atoms are remaining. We do not ob- 
serve heating or loss of molecules as induced by the purifica- 
tion sequence. 

We then use the pure, trapped sample of tunable halo dimers 
to study binary collisions. We measure inelastic decay, result- 
ing from the relaxation into more deeply bound states. In such 
a process the conversion of internal into kinetic energy by far 
exceeds the trap depth and leads to immediate trap loss of all 
particles involved. All our experiments are carried out in a 
regime of very low temperatures (k^T <C E\y), where the ini- 
tial kinetic energy of the colliding dimers is not sufficient to 
break up the molecules, as observed for 6 Li halo dimers in 
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FIG. 3: (color online) Scattering length dependence of the relaxation 
rate coefficient at 120 nK. The solid line is a linear fit to the data 
in the region a > 500ao (see text). The error bars refer to the statis- 
tical uncertainty. The shaded region indicates the a < r v( jw regime; 
here the experimental uncertainties are larger because the molecules 
have to be transferred through several avoided crossings with higher 
partial- wave states II 1511. 



Ref. (260. Moreover, spontaneous dissociation observed for 
85 Rb halo dimers |l4i [27|] is not possible as there is no en- 
ergetically open channel. Other density-independent losses, 
such as background collisions or light-induced losses, can 
also be neglected under our experimental conditions. We can 
therefore completely attribute the observed losses to inelastic 
dimer-dimer collisions. The decay of the trapped dimer sam- 
ple is thus described by the usual rate equation N = — oCj-qiRN. 
Here N indicates the number of dimers and a^i the relaxation 
rate coefficient. The mean molecular density n is given by 

n = [mcd 2 / (27tksT)] 3 ^ 2 N with m being the atomic mass and 
(Q denoting the geometric mean of the trap frequencies. 

We measure G^l as a function of the scattering length a for 
a fixed temperature T = 120 nK. We ramp the magnetic field 
to a desired value, and we then perform a lifetime measure- 
ment on the trapped dimers for storing times up to 2 s. As an 
example, the time evolution of the dimer number in the optical 
trap is shown in Fig. Ob) at two different values of the mag- 
netic field. We observe the expected non-exponential decay of 
the dimer number. We extract the value of G^ei by fitting the 
data with the above rate equation. The lifetime measurements 
are then repeated at different values of the magnetic field in a 
range from 7 G to 50 G. 

The observed dependence of the relaxation rate coefficient 
on the scattering length reveals an interesting behavior; see 
Fig. [3] Non-halo dimers exhibit a relative large and essentially 
constant collisional rate coefficient. In contrast, when the 
dimers enter the halo regime with increasing scattering length 
(a > r v dw), the rate coefficient first drops to a minimum. The 
minimum is found at a « 500<2o (~5r v dw)- F° r larger values 
of a, the rate coefficient increases with a. General consider- 



ations |28l I22I [30I] suggest to first approximation a universal 
linear scaling law for a^i according to a^i = C(h/m)a with 
a dimensionless constant C. From a linear fit for a > 500^o 
(solid line) we obtain an estimate for C of about 3. 

An essentially constant relaxation rate coefficient for vary- 
ing binding energy or magnetic field has been previously mea- 
sured in collisions between non-halo dimers. Class of the non- 
halo dimers includes 23 Na 2 [31], 87 Rb 2 I22D, 133 Cs 2 m various 
molecular states HHl^], and also 6 Li 2 /?-wave molecules as 
a process involving four identical fermions ll35ll . Our halo 
dimers composed of bosonic atoms thus show a novel and 
qualitatively different behavior. 

In a second set of experiments, we study the temperature 
dependence of the relaxation rate coefficient. In the ultracold 
domain, inelastic two-body collision processes are usually de- 
scribed in terms of a simple rate constant, i.e. a rate coef- 
ficient being independent of the particular collision energy. 



This applies to the case of inelastic atom-atom collisions [36], 



as well as to collisions between deeply bound dimers fl3_ 
In contrast to this usual behavior, we find a strong temperature 
dependence of the loss rate coefficient of halo dimers, so that 
a simple rate constant model does not apply. 

We have focused our measurements on the temperature- 
dependence in three different cases: the non-halo regime 
(120 ciq), the loss minimum in the halo regime (500 oq), and 
a more extreme halo case with increased loss (800 ao). As 
shown in Fig.lU the dimers exhibit the expected constant re- 
laxation rate outside of the halo regime (^9x 10- 11 cm 3 /s). 
In the halo regime, we observe a clear decrease of the re- 
laxation rate with decreasing temperature, both at 500 ao 
and 800 ao, roughly following a \[T -dependence. This sur- 
prising behavior raises the question whether the temperature- 
dependence of the relaxation rate is a property unique to halo 
states or whether it can also occur for other weakly bound 
Feshbach molecules. In the latter case, a halo dimer may just 
be seen as an extreme case of a weakly bound dimer. We spec- 
ulate that this observation is related to the fact that collisions 
between weakly bound dimers can involve more complex pro- 
cesses which go beyond simple two-body mechanisms. For 
instance, the release of binding energy may lead to a frag- 
mentation with three particles in the exit channel, i.e. a more 
deeply bound molecule and two free atoms. Breakup thresh- 
olds may manifest themselves in a more complicated depen- 
dence on the collision energy. 

The possibility to suppress inelastic loss by controlling tem- 
perature and magnetic field leads to a favorable situation, in 
which halo dimers exhibit a high degree of stability. In partic- 
ular, at 40 nK and 500 ao, we measure a loss rate coefficient as 
low as ~ 7 x 10 -12 cm 3 /s. This unusually small value corre- 
sponds to an order of magnitude improvement in the stability 
against collisional decay with respect to previously investi- 
gated cases of 23 Na 2 J3JJ] and 87 Rb 2 JH. 

To conclude, we have studied inelastic collisions between 
tunable halo dimers composed of bosonic atoms. We have 
observed a pronounced scattering length dependence, with a 
minimum in the loss coefficient as a most striking feature. We 
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FIG. 4: (color online) Temperature dependence of the relaxation rate 
coefficient a^i at 120 ao (triangles), 500 a® (squares), and 800 a® 
(circles). The solid lines are introduced as guides to the eye. 



have also found that inelastic loss is further suppressed by 
decreasing the temperature. The existence of the minimum 
raises the question whether this feature can be understood in 
terms of universal four-body physics, similar to a minimum in 
three-body recombination at large positive scattering length 
Q, SIH, H, which results from the destructive interference 
of two decay channels. The slow inelastic decay near the mini- 
mum may provide us with a favorable situation to study elastic 
dimer-dimer interactions or to search for universal four-body 
bound states |2D. 
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